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Abstract 

The relativistic equilibrium velocity distribution plays a key role in describing several high-energy 
and astrophysical effects. Recently, computer simulations favored Jiittner's as the relativistic gen- 
eralization of Maxwell's distribution for d = 1,2,3 spatial dimensions and pointed to an invariant 
temperature. In this work we argue an invariant temperature naturally follows from manifest 
covariance. We present a new derivation of the manifestly covariant Jiittner's distribution and 
Equipartition Theorem. The standard procedure to get the equilibrium distribution as a solution 
of the relativistic Boltzmann's equation is here adopted. However, contrary to previous analysis, 
we use cartesian coordinates in d+1 momentum space, with d spatial components. The use of 
the multiplication theorem of Bessel functions turns crucial to regain the known invariant form of 
Jiittner's distribution. Since equilibrium kinetic theory results should agree with thermodynamics 
in the comoving frame to the gas the covariant pseudo-norm of a vector entering the distribution 
can be identified with the reciprocal of temperature in such comoving frame. Then by combining 
the covariant statistical moments of Jiittner's distribution a novel form of the Equipartition The- 
orem is advanced which also accommodates the invariant comoving temperature and it contains, 
as a particular case, a previous not manifestly covariant form. 
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I. INTRODUCTION 



Incorporating the relativity principles in kinetic theory is crucial not only to understand 
the theoretical grounds in the description of relativistic many particles systems [IJ [2] but to 
interpret relativistic high-energy experiments like those involving heavy-ion collisions j3] as 
well as phenomena in the astrophysical jl] and cosmological realms [5]. These include, for 
instance, the use of the relativistic Bolztmann equation to understand the thermal history 
of the Universe [6l [7] and the structure of the Cosmic Microwave Background Radiation 
spectrum associated to its interaction with hot electrons in galaxy clusters [8]. 

In the case of equilibrium the history of the relativistic analogue of Maxwell's velocity 
distribution goes back to F. Jtittner, who in 1911 turned to relativity to consistently get rid 
of the contribution of particles with speeds exceeding that of light in vacuum, denoted by 
c, and which are contained in Maxwell's distribution. This was achieved upon replacement 
of the non-relativistic energy of the free particles by its relativistic form and a maximum 
entropy principle thus yielding the so called one-particle Jiittner's distribution [U] 

n Vp2 c 2+m2c4 

J A-nkTm 2 cK 2 {mc 2 /kT) ' 1 ; 

which is written here in momentum rather than velocity space, m is the rest mass of 

any of the particles forming the gas, k is Boltzmann's constant and K 2 is the modified 

Bessel function of order two |10j . The remaining quantities: particle number density n, 

temperature T and the relativistic spatial momentum of the particle p, should be here 

understood as determined by an observer comoving with the gas as a whole. The mass shell 

condition, (E/c) 2 — p 2 = m 2 c 2 , is assumed along the sequel. In the non-relativistic regime 

|p| -C mc,mc 2 <C kT, we have / ~ /Maxwell, where /Maxwell stands for Maxwell's velocity 

distribution. Thus particles's speeds beyond c are just an artifact of the non-relativistic 

approximation. 

Jiittner's distribution in the form can be disadvantageous in that it is not manifestly 
covariant, namely it is not expressed in terms of Lorentz tensors which in turn explicitly 
show its behavior under Lorentz transformations to investigate its description for frames in 
relative motion. To make it manifestly covariant two key information pieces are needed: the 
transformation under change of frame of / and that of T. Both of them have been considered 
in the literature at large [Hi [12] and [TJ [131 US US EL f° r temperature and distribution 
function, respectively. 
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Upon multiplying Q by the characteristic function for the box confining the gas so that 
X Box (x) = 1 for x within the box, and zero otherwise, the resulting distribution /(x, p) := 
X Box (x)/(p) is defined on the single-particle phase space. An observer would thus determine 
/ particles of the gas in a volume d 3 x located at x and having momentum p within range 
d 3 p. Moreover / must be Lorentz invariant [131 EE Ej • This is readily seen in the case of 
equilibrium [18J: for a simple gas the number of particles M is invariant and so jV = J dfif, 
with dfi = d 3 xd 3 p, must be. Since dfi is a Lorentz invariant measure due to a compensation 
between the transformations of the (mass-shell) spatial momentum and space measure, / 
must be so too. Now since x Box * s invariant [44] then / is. Hence a manifestly covariant 
form of / amounts to manifest invariance and, in particular, it should involve the behavior of 
temperature under Lorentz transformations. The manifestly invariant Jiittner distribution 
was given long ago [21 [HI [151 US]- It was determined by introducing spherical coordinates 
on the pseudosphere associated to the mass shell in relativistic four-momentum space. This 
provided a rather elegant and fruitful treatment which however does not seem to have been 
fully appreciated in particular regarding temperature. 

Alternatives to ([I]) have been proposed recently for which Lorentz covariance is incorpo- 
rated in a different manner. In [20] a quantum mechanical approach including an invariant 
"historical time" was considered to derive a manifestly covariant Boltzmann equation with 
stationary solution implying a different ultra-relativistic mean energy-temperature relation 
whereas in [2H .22] a maximum relative entropy principle was combined with an invariant 
group theoretical measure approach to obtain an equilibrium distribution. Both alternatives 
have been recently criticized in [23]. Moreover recent Molecular Dynamics Simulations of a 
two component one dimensional simple relativistic gas showed agreement with ^ and tem- 
perature defined through an equipartition theorem was shown to be invariant under change 
of frame [23]. The study of the two-dimensional case along these lines has been reported in 
[25] and [2H]- For three spatial dimensions Monte Carlo simulations have been considered 
favoring also Jiittner's distribution [2?]. Amusingly, as for temperature, such kind of analy- 
sis take us back to the long standing discussion of whether a moving object appears cooler 

PHE21 122 ESI . 

In this work we shall obtain the manifestly invariant Jiittner distribution by adopting 
"cartesian" rather than spherical coordinates in relativistic (rf+l)-momentum space [HI 
[T5] . Our approach can be considered as an alternative to such previous derivations of the 
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Jiittner's distribution and to others [TJ |28] in the sense that we avoid adopting any frame 
along the sequel which otherwise would raise the question on the Lorentz transformation of 
temperature. We only allude to the frame comoving with the gas at the end of the analysis 
to relate the kinetic description with thermodynamics in particular to identify temperature 
with the Lorentz invariant norm of a thermal vector. We once more obtain in this way 
the thermal four vector introduced long ago [21 [2J [EJ [19] which is formed by the product 
of the inverse comoving temperature with the four velocity of the g whole. Also 

the lower dimensional cases recently studied [2H ESI ESI EZ] are contained in our results. 
Hence comoving temperature is seen to play a role analogue to rest mass of a particle. The 
compatibility of such an interpretation is further investigated in relation with a manifestly 
covariant form of the equipartition theorem. 

The paper is organized as follows. For the sake of completeness section II briefly re- 
views the derivation of the Jiittner distribution as an equilibrium solution of the relativistic 
Boltzmann equation. This includes a brief description of the two types of approaches: the 
manifestly covariant one adopting spherical coordinates on the mass shell pseudo-sphere 
in momentum space and the one adopting the comoving frame and cartesian components. 
Section III provides the details of our analysis in which we use "cartesian" coordinates in 
momentum space to get the manifestly invariant distribution. In particular the thermal 
Lorentz vector is here characterized. Its role in the relativistic covariant equipartition theo- 
rem is the subject of section IV. Finally section V summarizes our results. We also describe 
the behavior of a Planckian distribution when use is made of invariant temperature. The dif- 
ficulties to define a non-comoving temperature for a gas of massive particles is also touched 
upon to further stress that invariant comoving temperature seems to be the only consis- 
tent possibility to define temperature according to the standard relativistic kinetic theory 
framework. 

Our conventions are the following. Spacetime components of tensors are denoted by greek 
indices: fi, u, ■ ■ ■ =0,1, 2, 3, the zero component referring to time whereas spatial components 
are denoted by latin indices i,j,--- = 1, 2, 3. Einstein sum convention is assumed for both 
types of indices throughout. The Minkowski metric has components rj^ = diag{+, —,—,—}. 
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II. JUTTNER DISTRIBUTION AND THE RELATIVISTIC BOLTZMANN 
EQUATION 



Consider a simple relativistic gas composed by identical particles of mass m. To each 
particle correspond space-time and four-momentum coordinates, respectively given by 
and p^. The evolution of the one-particle distribution function of the gas is governed by the 
relativistic Boltzmann equation [TJ |28l [29] 



|7 = 



C[f*J] := / (./v./--/,/) Fadn^. (2) 
J Pi 



Here C[f*,f] is the so called collision term and /* = /(x, p*,t), * implying a quantity 
is evaluated after the collision, cr is the differential cross section of the scattering process 
whereas Q is the solid angle. /C M denotes an external 4-force while F is the so called invariant 
flux, F = a/ (PiPh) 2 — m^c 4 , and ^fr the invariant measure on mass shell. We shall be 
interested in the case /C M = 0. 

The equilibrium state can be defined so that the entropy production vanishes [U [HI [15J, 
[28] . In covariant form this means = 0, with the entropy flux given by 

S» = -kc J ^f\n(h 3 f), (3) 

and h is a constant with dimensions to make the argument of the logarithm dimensionless. 
Zero entropy production requires f*f*i = ff\ and so the collision term C[f*,f] becomes 
zero too. Such condition can be written as 

lnr + ln/ 1 *=ln/ + ln/ 1 , (4) 

where x Box canceled out. Eq. Q is fulfilled by the collisional invariants which for the 
present case become the four-momentum up to a constant (independent of so that an 
/ solving (|4]) takes on the form 

In / = - (A(x) + 9,(1)/) & f = a(x) exp (-6 M (x)p^). (5) 

Here the negative sign has been introduced for later convenience, is the 4-momentum of 
a particle of the gas and a = e~ A . The equilibrium distribution function is fully obtained 
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by requiring consistency between ^ and the left hand side of ^ equated to zero. Such 
substitution yields A independent of x M and d u Q^(x) + d^Q v {x) = 0, whose solution is 
Qfj,(x) = oj^ v x u + 6 M and uj^ v = —u Vfi . These correspond to the ten Killing vectors (6 
for uo^ v x u and 4 for B M ) under which Minkowski spacetime metric is invariant. They are 
associated with Lorentz transformations and spacetime translations, respectively [30J. As 
we will see below B M can be identified with the four velocity of the fluid as a whole and thus 
it inherits spacetime symmetries in the form of rigid motions: world lines of neighbor fluid 
elements would keep their separation whenever they lie along Killing vectors |14j . We shall 
restrict such motions to translations in the present work and so only M is considered. 

To determine a and G M one assumes that typical physical quantities are related to the 
statistical moments of the distribution. For instance the particle number density flux and 
the energy-momentum tensor 

N» = cJ ^VV , (6) 



c I PV7^ , (7) 



can be respectively written as 



= ~ac§- , (8) 



from which the denomination of I as generating functional suggests itself [14] . Having 
explicitly defined I will allow to determine a. On the other hand to obtain M one relates 
the kinetic theory form for thermodynamical quantities with equilibrium thermodynamic 
equations assumed to hold in a frame comoving with the gas as a whole. This will relate G M 
with a comoving temperature. 



To evaluate ( 10 ) one can express it in spherical coordinates on the mass shell in momen- 
tum space [TU HS1 [12], namely, since p° = a/p 2 + m 2 c 2 , one has for the components of the 
particle's momentum 

= (jnc cosh x, mc sinh x sin ^ cos 0, 

mc sinh x sin ^ sin 0, mcsinhxcos#) , (11) 

0<X<OO, 0<6 <7l, < </> < 2tT , 



on the pseudosphere = m 2 c 2 . Hence (10) becomes 



J = / (ift (3) (ixsmh 2 xe- mcecoshx 



where dil^ is the element of solid angle in three spatial dimensions, K\ is the modified 
Bessel function or order one pi)] and M has been chosen to lie along the x — axis. Upon 
use of the identity between Bessel functions 4-[u~ n K n (u)] = —u~ n K n+ i(u) to evaluate the 
components of the statistical moments ^ and (|9| allows to relate their components in the 
comoving frame to get the equation of state for the gas V = pc 2 0, with V the pressure and 
p the density of energy, both in the comoving frame. This led Israel [2] to identify = -§j,, 
with T the invariant comoving temperature. 



A different way to deal with (10) is to consider that, being invariant, it can be calculated 
in a convenient frame, say one in which B^ 1 = (0°, = 0). Here cartesian coordinates have 
been used in momentum space [H [28] . Noticing that (i) M in (|6| should be timelike to make 
the integral to converge and (ii) the only available timelike vector for the gas as a whole is 
its velocity U^, one is led to propose M = so that T is identified with a quantity in a 
frame in which IX, the spatial part of IP, is zero. Such frame is the comoving frame to the 
gas. This approach in which one picks the comoving frame from scratch begs the question 
about which is the Lorentz transformation of the temperature. 

It would be desirable to be able to combine the power of the spherical components of the 
manifestly covariant approach mentioned afore with the more intuitive and easier to handle 
cartesian ones and still be able to investigate the behavior of temperature under Lorentz 
transformations. This is indeed possible and will be discussed in the next section. 



III. MANIFESTLY COVARIANT GENERATING FUNCTIONAL WITH CARTE- 
SIAN COORDINATES 

Since our analysis remains the same independently of the number of spatial dimensions 
we consider such general case at once. Hence unless otherwise stated, tensor indices run as 
follows: p:, v ■ ■ ■ — 0, 1, 2, . . . , d, for d spatial dimensions. 
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A. Determination of a 



Let us consider a frame non-comoving with the gas so that the spatial components of 



G M , 0, are non-zero. In <i-dimensional space the integral X, Eq. (10), just requires to 
change the measure from d 3 p to d d p. Now we adopt spherical coordinates only for the spa- 
tial part of the momentum and assume that © • p = |0||p| cos$i. We then have d d p = 
Ipl^dlpltZfiW, where dflW = (sintfi) d - 2 (smtf 2 ) d ~ 3 • • • (sintf^) 1 ^!^ • • • d-& d ^ 2 d(p = 
Yl d Zi(s™$ i ) d - i - 1 d'& i d<p, where 0<<p<27ry0<tfj < 7r [31 j . Thus X can be written 
as 

X=Sd _ x j rf IPllPl d " 1 (sini?i) d-2^ ie -e p° e |e||p|costf lj (13) 



d-1 



Here S^-i = 2 l d t x \ is the hyper-surface of the d — 1 dimensional sphere [31] resulting from 
integrating over and dQ^' 1 ^, which excludes To integrate over $i it is better to use 
the series form of the spatial exponential e' e " p ' cos1?1 = Y2T=o ^ e ^ p ^ c ! ° stfl ' > . In this way (13) 
becomes 

\@\ k f |p| d - 1+fc d|p| 



] n ** 



k=0 

X 



■_ e -e oP o x 



P7T 

/ sintf^costf^cMi. (14) 
Jo 

The remaining angular integral is non-zero only for k = 2n, n = 0, 1, 2, . . . and, moreover, 
it can be related to the beta functions B (^-^, [10J so that 

2n+d-2 

D °[p° 2 -m 2 c 2 ] 2 dp . (15) 



x / e 



Where we have exchanged the independent variable |p| by p° by using the mass shell con- 
dition. The further change y° = p°/mc and Zq = tucQq together with the integral form 
dxe~ ax (x 2 — l) m ~5 = (|) m T ^™t^ K m (a) for the modified Bessel function [10] leads to 



T - <? vl©| 2 " ^2n + l d-l\ 



' mc )2n+d-l x 



n=0 

2n+d-l 



x(-J ^m^ ^^(^o). (16) 
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At this point we make the following considerations: (i) We use the known expression of 
the beta function in terms of the gamma functions JTU] and then (ii) introduce z^ = mcQ^ 
together with (3 := \z\/z°, < (5 < 1, the latter following from the fact that M is timelike. 



Thereby Eq. ( 16 ) may be rewritten as 



d+l 



X = 2 2 (mc 



\d-l 



71 



°° K„,d-i(zo) 

n=0 



2 n n\ 



(17) 



Here we arrive at a critical point from the technical perspective. The sum in (17) can be 



further reduced to a modified Bessel function upon using the multiplication theorem of the 
Bessel functions 132 



m=0 



with I A 2 — 1| < 1. This finally produces 



1= 2(mc) d - 1 2 



ml 



d-l 
7T\ — 



(18) 



K^i (z) 

2 



(19) 



Formula ( 18 ) is essential to trade the Bessel series in ( 17) by a single modified Bessel function 



in (19). This becomes (12) when d = 3. We see that X is only function of the invariant z. 



Now we can obtain the relativistic particle number density flux from ([8l producing 



N» = 2m d c d+1 a ( 2 



Such equation can be solved for a giving 



7T 



N 



a 



Kd+i (z) — 

2 Z 



d-l 

/ mcO\ 2 



(20) 



2c(mc) d K d±1 (mc6) \ 2n J 



(21) 



It is clear from (20) that z^ and iV M point in the same direction and since = NU^/c we 
have that 

ft 

(22) 



9" = —W*. 

c 



The physical significance of within our approach is discussed in the next subsection. 



B. 0^ and the invariant Temperature 



Now we seek to identify M with a thermodynamical quantity. Let us take as our starting 
point the Gibbs form of the second law of thermodynamics for a closed system [151 
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which we shall assume to hold in the comoving frame 

SU = T5S - V5V. (23) 

Clearly we must relate the relativistic kinetic expressions like the energy-momentum tensor 
and entropy flux with internal energy, U, entropy, S, pressure V and volume V appearing 
in (23). To begin with let us introduce the distribution function ^ and Eq. (21), in the 



expression for the energy-momentum tensor (J9|. This leads to 

^--e^-i^-e-™)- (24) 

The comoving pressure can be obtained from the energy-momentum tensor in ^-dimensions 

as 

1 N 

v = ~d A ^ T " u = e ' (25) 

where the projector is = 7]^ u — IPLPc -2 . The corresponding <i-dimensional entropy flow 
(|3| can be conveniently reexpressed as 

= -k [In (ah d )N^ - T^Q U ] . (26) 

It is worth stressing that the quantities iW, T^ u and are densities, and therefore do 
not depend on the size of the system. It is however rewarding to include the fact that the 
gas is inside a box and to do so we make use of the characteristic function x Box defined in 
the introduction. In particular integrating over d spatial dimensions 



•A/" = c 1 J X^N»da„ (27) 
S = c- 1 j x^S'da,. (28) 
= c' 1 j X Bo T^d*„, (29) 



Q 



In the previous expressions (27)-(29), do^ = n^d I S, with a unit normal to the spatial 
hypersurface, and therefore a timelike vector; in particular d d S = d d x when = 5q. 
Obviously TV and S are Lorentz invariants while is a vector. The latter is just the total 
relativistic momentum of system. 



Combining (26) and (28) yields 

S = -k [AT In (ah d ) - 6^] , (30) 
11 



Now we can make contact with (23) by going to the comoving frame. Consider (30) for 



which such process means, in light of (22), M £ M 



QQ°. The differential of the 



entropy, (30), in thermodynamic space (G constant) becomes hereby 



5S = kO 

kQ 

c 



5Q U - 

c N\ 

[5(cG°) + V5V] 



(31) 



where use has been made of (21), (25) and the relation 5V 



AfSN 
N 2 ■ 



with (31) gives rise to the identification 



6 



c 

kf' 



Comparison of (23) 



(32) 



with T the comoving temperature by interpreting cQ° as the relativistic generalization of 
internal energy U; this clearly holds in the low speed regime in which cQ° ~ A/"mc 2 + f/ non _ re i, 



with ?7 non -- re i the usual non-relativistic internal energy of the gas. Thus (21), (22) and (32) 
complete our analysis of the manifestly covariant distribution function which is expressed in 
terms of invariant quantities, N, 0, m, c, O^p^, in the fashion 



N 



(^) ' «pW)- 



(33) 



2c(mc) d K d+i (ttccO) 

Eq. Q is obtained from (33) by considering d = 3 and using the comoving frame to the gas. 
It should be stressed however that in deriving (33) no assumption was made in regard to the 
Lorentz transformation of temperature. The latter came from the requirement to reproduce 
equilibrium thermodynamics in the comoving frame. 



IV. COVARIANT EQUIPARTITION THEOREM 

Equipartition theorems are usually considered to connect temperature with averages of 
say kinetic energy in the non-relativistic case and certain functions of momenta in the 
relativistic case [121 [MJ. Now, as we have shown in the previous section a manifestly covariant 
approach to relativistic equilibrium distribution function unveils the convenience of using an 
invariant temperature which is the one associated with the comoving frame of the gas. Thus 
it is of interest to investigate the role an invariant temperature plays within a manifestly 
covariant equipartition theorem. 
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The non manifestly covariant but relativistic equipartition theorem seems to have been 
first considered by Tolman [31] and later by others |35J . Its use as a criterion to determine the 
Lorentz transformation of temperature was criticized by Landsberg [T2] who stressed that it 
can accommodate both an invariant as well as a transforming temperature under a change 
of frame. Recently Cubero et al. performed numerical simulations indicating the existence 
of an invariant temperature [21] on the basis of the relativistic equipartition theorem of [12] 
thus hinting at the invariant temperature included in Landsberg's analysis. In this section 
we shall provide a manifestly covariant form of the equipartition theorem which not only 
contains an invariant temperature but includes Landsberg's version. Moreover our form of 
the theorem is neatly expressed in terms of the total momentum of the gas. 



Let us start by reexpressing the generating functional X (10) in (d + 1) -dimensional 
momentum space [1] 



(34) 



where H(p°) is the Heaviside function that restricts (34) to positive energies. The mass shell 
condition is accounted for by the Dirac's delta. 



By covariantly extending the argument in [34] , we next integrate by parts d + 1 times 
-one for each p^- and discarding the boundary terms at infinity we obtain 

2 



1 



d + 



~ [p u w- \e' & ^ S (p a p° - m 2 c 2 )] H(p°) d d+1 p . 
1 J op v 



(35) 



To integrate over p° we apply the properties of the Dirac delta which finally yields 



IP] 



d d p 
"po"' 



(36) 



Remarkably, notwithstanding the second term in (36) containing the spatial components pt 
it is Lorentz invariant overall. 



To connect with the physical framework Eq. (36) is plugged back into (j8|) and ^ to 
yield 

2 



N c 



c f dp 



d J p° 
c f d d p 



f<P° 



IPI 
p° 



dp i 



■ dp 



dp 1 



Ipj 
p° 



dp* 



dpi 



-p' 



(37) 
(38) 



13 



A further spatial integration of (37) and (38) according to (27) and (29) gives rise to 

ftp 1 * ' 



d = e 



d 



dp 1 



e, 



op 1 



: dp a 
dp 1 



(39) 
(40) 



where ((•)) := j? j ■ d d \i has been used. It is worth stressing that Eq. (39) is just the 



manifestly covariant form of the equipartition theorem corresponding to that advanced by 
Tolman and Landsberg [121 EH] , respectively, expressed using the invariant comoving tem- 



perature T = c/kQ, Eq. (32). A comment regarding covariance of (39) is here in order. 



Although it contains the spatial sum p % -§^-, we arrive to it from the manifestly covariant 



equation (|34|). This is analogous to recognizing the invariance of Clearly the argument 



holds also for both terms in the r.h.s. of Eq. (40). 



Observing that (39) and (40) contain a common term suggests that the relativistic mo- 



mentum can be made to enter the covariant equipartition theorem. To proceed further we 



need to determine the first term in (40). First we combine (29) and (24) to obtain the 
explicit form for the momentum 



Q a =N 



62 



Kd+s (mc6) 
mcQ— - — ; — — 1 



Ka 



(mcG) 



(41) 



By equating the projections Q a G a of (|40|) and (|41|) we identify 



d 



p 



■ dp u 

dp 1 



V 



Kd+3 (mcB) 

1 mcQ. 



From here we finally arrive at 



— z 



r d {z) 



Kd+i (mcQ) 



J-'d(z), z = mcQ 
Kd+3 (z) 

-1-z 



(42) 

(43) 
(44) 



Kd+i (z) 

2 

Amusingly just in the same way as the energy of a particle having momentum Pp article is 
determined by an observer having velocity U% hs is given by E obs = Pi/ part icie^obs we interpret 
as the energy of the gas determined by the comoving observer. 
Td{z) reduces to the usual values of the equipartition theorem in the limiting cases (See 
Fig. §. 

- z 1 

TA,) < 2 ' >>i; (45) 
d, z < 1. 
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Figure 1: The graph shows J~d(z) vs. z = mc 2 jkT for d = 3. J~s{z) corresponds to the quotient of 
the average of the relativistic kinetic energy per particle and the comoving temperature. We see 
that for z > 1 and z < 1 we get the appropriate limits, while for intermediate values of z, the 
ratio is a fraction between 3 < ^(z) < 3/2. 



Notice that although a relativistic equipartition theorem had been considered before by 
other authors [I2J ES] their approach was not manifestly covariant. Their work and ours 



coincide in terms of the function Td (44), in particular the behavior in the non-relativistic as 
well as in the ultra-relativistic limits. Remarkably, recent numerical calculations adopting 
Monte Carlo methods [27] confirm Td as giving the relativistic kinetic energy divided by kT. 



V. DISCUSSION 



The interest in incorporating the relativity principles into kinetic theory goes beyond the 
theoretical foundations: Actual observations and experiments like for instance in high-energy 
physics [3J, astrophysics jl] and cosmology [3] require a description of relativistic many- 
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particles systems in terms of, say, Boltzmann's equation and the corresponding equilibrium 
distribution [HI HE]- Recently, Cubero et al. [24J (See also references there) have developed 
numerical simulations based on Molecular Dynamics pointing to Jiittner's as the equilibrium 
distribution in agreement with their numerical analysis, as opposed to other proposals in 
the literature, for one and possibly other number of spatial dimensions. This was further 
confirmed for two j25j [26] and three [27] spatial dimensions, the latter adopting Monte 
Carlo simulations instead. In [23] , the old problem regarding the relativistic transformation 
of temperature [UJ was also considered in connection with a previous relativistic version of 
the Equipartition Theorem proposed by Landsberg [T2]. Based upon such theorem Cubero 
et al. determined, within their approach, that temperature should possess an invariant 
character. It is worth mentioning that such analysis were not formulated in a manifestly 
covariant form and indeed some further insight was needed to interpret the contributions 
entering such theorem as well as temperature. 

In this work we have obtained a new derivation of the manifestly invariant Jiittner's 



relativistic distribution function (33). This is based on cartesian coordinates in 
momentum space (d spatial dimensions) in contrast with the known results developed using 
spherical coordinates. This was made possible by the use of the multiplication theorem for 



Bessel's functions ( 18 ) which simplified the treatment of a series involving Bessel functions 
|32j . In this approach no assumption is made a priori of any specific relativistic character for 



temperature. The latter appears through the invariant norm, G, of a four- vector, M (22), 
and it is invariant just for the same reason a point particle's rest mass is. Indeed because 
of the assumption that relativistic equilibrium kinetic theory should agree with standard 
thermodynamics for an observer comoving with the system under study the pseudo-norm 



becomes G = c/kT, Eq. (32), T being the comoving temperature of the gas. Finally we have 



developed a manifestly covariant Equipartition Theorem, Eq. (43), in which the average of 
the energy-momentum of the gas as determined by the comoving observer is given by a 



function jF d , Eq. (44) of the invariant temperature. Indeed this is analogous to the case 
of a point particle for which the energy is obtained by projecting momentum along the 
four-velocity of the observer. Here we have the thermal vector M = with K/ 1 the 

four-velocity of the gas as whole which defines a comoving observer reading the invariant 
temperature T. A further comment is here in order regarding the difference between our 
approach and previous ones. While previous versions of the Equipartition theorem 



16 



relate temperature with a peculiar combination of relativistic quantities (See for example 



Eq.(39)), here we interpret that combination as included in (43) which relates the invariant 
temperature with the averaged relativistic momentum, Q 1 * . 

In a nutshell the manifestly covariant form of Jiittner's distribution leads naturally to 
consider the invariant comoving temperature to characterize the equilibrium regime. As an 
illustrative example of how this result applies to known cases let us consider the case of black 
body radiation. Let us recall the standard analysis of this problem [5j [7J [18] : The Lorentz 
invariance of j^, containing the specific intensity I v and the frequency v of the photons, 
implies the invariance of the Planckian distribution 

1 (46) 



2hv 3 evr-i 
Indeed for two frames in relative motion one should have 

hv hv' 



(47) 



kT kT r 

where primed quantities refers to an observer that moves with respect to radiation and, in 
particular, T 7 is suggested as the non-comoving temperature. Since the photons suffer a 
Doppler shift: 

^ = 7(1-/3 cos 5), (48) 

where 5 is the angle between the momentum of the photon and the gas velocity whereas (3 
is the ratio between the gas speed and c. The quantity T takes the anisotropic form 

T 

7(l-/3cos«J)' 1 } 



Notice that we could alternatively write a manifestly invariant form for (46) as (e — l) -1 . 
In this way one can focus in the invariant product O^jA By evaluating it in the two frames 
mentioned above and using that for photons p'° = |p'| 

^ = 9'V°(l-/5cos5). (50) 

Considering the invariant comoving temperature T the components 0' M become ^ (c, U) and 



thus ( 50 ) reduces to ( 48 ) . This shows the consistency of adopting such invariant comoving 



temperature, trough M , in the description of the black body radiation without resorting to 



T', Eq. (49). The latter can be regarded only as an auxiliary quantity for the following rea- 



sons. Firstly, actual observations of the Cosmic Microwave Background Radiation (CMBR) 
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reveal there is a frame in which it presents the black body structure. Measurements 
however involve brightness depending on frequency rather than a non-comoving temperature 



(49). Secondly, in the case of massive particles one finds and obstacle to handle (49) (See 
for instance [38J). For massive particles we would have 

T 

T = 7 (1 -(3 p f3 cos 5) ^ 

a = M 

which makes no sense due to the dependence on the particles's momentum. However adopt- 
ing an invariant comoving temperature is viable just for the same reason that the case for 
photons described above works. 

Our work then adds to recent claims elaborated on the basis of an Unruh-DeWitt detector 
[39] pointing to the impossibility to have a relativistic transformation of temperature [40J. 
In any case all these results reinforce the idea that temperature makes undisputable sense 
in the comoving frame. 

There are several possible extensions of the present work which could be of interest. One 
possibility is to extend the analysis in the present work to the case of relativistic Fokker- 
Planck equation, say along the lines of [JT]. In other direction, we have seen that in order 
to solve Boltzmann's equation for the equilibrium distribution the concept of rigid motions 
appeared. They can be seen related to the symmetries of Minkowski spacetime that we have 
here considered and specifically to the existence of Killing vectors. Although symmetric 
curved spacetimes have been studied previously for instance in j2j [6] there seems to be no 
attempt to incorporate the analogue of Noether's theorem in regard to the kinetic theory in 
particular noticing that the latter can be linked to relativistic quantum field theory pfl E3] . 
This would be particularly useful in the manifestly covariant approach. Generally covariant 
approaches to statistical mechanics have also been studied to investigate many particle 
systems in general relativistic theories (42] • It is a challenge to incorporate similar ideas to 
kinetic theory. Furthermore the very notion of spacetime has been considered in these terms 
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